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Abstract 

This  paper  presents  an  optimization  approach  applied  to  a  whole  fuel  cell  (FC)  air  supply  system  including  its  geometry  and  its  control.  The  aim 
is  to  optimize  its  power  consumption  along  with  its  mass.  Particle  swarm  optimization  (PSO)  algorithm  is  used  to  define  the  design  parameters 
of  both  permanent  magnet  synchronous  motor  (PMSM)  and  a  fuzzy  logic  controller  (FLC).  The  results  are  compared  with  those  obtained  by 
a  sequential  optimization  process  and  advantages  of  co-design  optimization  approach  are  clearly  shown.  Indeed,  a  significant  reduction  of  the 
objective  function  (made  up  on  both  motor  mass  and  energy  consumption)  on  a  considered  operating  cycle  can  be  obtained. 

©  2008  Elsevier  B.Y.  All  rights  reserved. 

Keywords:  PSO;  Fuzzy  logic  controller;  Permanent  magnet  synchronous  motor;  Co-design;  Fuel  cell  air  circuit 


1.  Introduction 

Since  90s,  problems  of  the  greenhouse  gas  emissions  con¬ 
trol  and  hydrocarbon  energy  resources  exhaustion,  started  again 
researches  on  fuel  cell  (FC)  in  many  fields,  such  as  transport,  sta¬ 
tionary  power  generation  and  portable  applications  [1].  Among 
these  fields,  the  surface  transports,  that  strongly  use  oil-based 
hydrocarbons,  is  one  of  the  sectors  generating  the  strongest 
greenhouse  gas  emissions.  On  the  opposite,  a  FC  directly  sup¬ 
plied  by  hydrogen  does  not  generate  locally  any  environmental 
pollution.  From  this  point  of  view,  it  is  an  interesting  alternative 
to  the  internal  combustion  engine. 

A  proton  exchange  membrane  FC,  which  is  the  most  com¬ 
monly  used  FC  in  the  transport  area,  is  an  electrochemical 
converter  which  delivers  electrical  power  and  thermal  energy 
by  a  redox  reaction  starting  from  a  fuel  (generally  hydrogen). 
The  energetic  efficiency  of  FC  itself  is  relatively  high  com¬ 
pared  to  other  more  conventional  technologies  (#60%  for  proton 
exchange  membrane  fuel  cell  efficiency  and  #32%  for  internal 
combustion  engine  efficiency).  But  its  output  power  decreases 
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due  to  ancillary  circuits  (fuel  (hydrogen)  and  fuel  oxidizer  (oxy¬ 
gen)  supply,  humidification,  cooling,  electric  converter,. . .)  that 
are  necessary  to  the  correct  operation  of  the  FC  (Fig.  1).  Among 
these  ancillaries,  the  air  supply  circuit  (oxygen)  of  the  FC  gen¬ 
erator  is  particularly  “greedy”  in  energy  [2].  It  is  classically 
carried  out  by  a  motor-compressor.  Electrical  conversion  net 
efficiency  of  the  energy  brought  by  the  fuel  is  around  50%  but 
up  to  35%  of  the  produced  electrical  energy  is  consumed  by 
the  ancillaries  (among  this  part,  the  2/3  are  consumed  by  the  air 
supply).  Thus,  air  supply  system  optimization  is  an  important 
milestone  on  the  road  of  efficient  FC  systems  [3].  In  this  context, 
the  main  goal  of  this  paper  is  to  propose  a  new  methodologi¬ 
cal  optimization  approach.  In  manners  of  electrical  engineering, 
the  classical  approach  is  to  use  sequential  optimization.  In  this 
case,  the  design  process  has  several  steps  starting  from  system 
geometry  optimization  to  finish  with  the  calculation  of  the  con¬ 
trol  laws.  Such  a  methodology  has  already  been  proposed  in 
several  papers.  In  [4]  the  driving  engine  design  of  the  com¬ 
pressor  unit  is  proposed.  Similar  works  are  carried  out  on  other 
engines  in  [5,6].  [2]  deals  with  an  optimization  of  the  air  supply 
circuit  controller  to  improve  system  performances.  More  gener¬ 
ally,  controller  optimization  for  many  engines  are  carried  out  in 
[7,8].  In  [9],  actuator  control  optimization  is  performed  based 
only  on  the  optimization  of  the  scaling  factors  of  a  FLC.  In  [10], 
PSO  is  used  to  define  parameters  of  a  PID  controller  in  order 
to  improve  its  performances.  But  is  it  really  the  optimal  system 


122 


A.  Sari  et  al.  /  Journal  of  Power  Sources  179  (2008)  121-131 


Hydrogen  tank 

r=~^- - •, 

H2  circuit  i 


Electrical 
current  if 


DC/DC 

converter 


tank 


Compressor 
air  circuit 


-4 


Humidififf 


UFC  stack 


Cooling  "1 
Circuit  i 

>  TO1 


Motor 


Turbine 


power  and  it  became  a  very  common  choice  in  driving  technol¬ 
ogy  over  the  last  few  years  due  to  some  of  its  inherent  advantages. 
Those  include  high  torque  to  current  ratio,  large  power  to  weight 
ratio,  high  efficiency  and  robustness. 

PMSM  is  modelled  with  an  equivalent  Park  machine  which 
makes  it  possible  to  simplify  motor  equations.  Inputs  of  this 
model  are  the  three  phase  voltages  (ua,  r>b>  fc)  coming  from  the 
inverter.  Thanks  to  Park  transformation,  forward  and  quadratic 
voltages  (ud,  uq)  can  be  obtained  from  these  three  voltages. 
Motor  flux  expressions  on  forward  and  quadratic  axes  make  it 
possible  to  evaluate  forward  and  quadratic  currents  (id,  iq)- 


Fig.  1.  Basic  schema  of  a  fuel  cell  system  [4]. 
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that  is  obtained  with  all  these  sequential  methodologies?  Our 
proposal  is  to  compare  this  classical  approach  to  a  simultaneous 
optimization  (i.e.  optimization  of  the  whole  system  at  the  same 
time)  methodology.  Such  an  approach  has  been  presented  in 
[11,12]  concerning  mechanics,  electronics  and  micromachine 
domain  and  is  similar  to  a  multi-objective  optimization  [13]. 
The  co-design  being  an  emergent  idea  in  electrical  engineering, 
the  references  are  still  limited.  Our  aim  will  thus  be  to  apply  a 
co-design  approach  in  order  to  evaluate  its  effectiveness  and  to 
compare  it  with  the  sequential  one. 

The  paper  is  organized  as  follows.  In  Section  2,  an  overview 
of  the  driving  motor  and  compression  head  modelling  is  pre¬ 
sented  and  fuzzy  logic  controller  is  introduced.  Section  3 
presents  particle  swarm  optimization  (PSO)  and  its  imple¬ 
mentation  on  the  considered  co-design  problem.  The  results 
obtained  by  both  optimization  methodologies  (sequential  and 
simultaneous)  are  compared  in  Section  4.  Finally,  Section  5  con¬ 
cludes  this  paper  by  pointing  the  key  results  and  the  possible 
outlooks. 

2.  Air  supply  circuit  modelling  and  control 

The  considered  air  supply  subsystem  is  here  a  motor- 
compressor  built  up  from  a  rotary  compressor  and  from  a 
permanent  magnet  synchronous  motor,  which  enables,  asso¬ 
ciated  to  efficient  control  laws,  to  supply  the  FC  stack  with 
the  air  flow,  in  order  to  produce  the  desired  electrical  power 
[14].  In  case  of  use  with  a  5kW  proton  exchange  membrane 
fuel  cell,  considered  specifications  for  this  system  are  the 
following: 

•  Maximum  mass  flow:  10gs-1, 

•  Maximum  absolute  outlet  pressure:  2  bars, 

•  Rated  rotational  speed:  3000  rpm, 

•  Nominal  mass  flow:  6  g  s_1 . 

2.7.  Analytical  model  of  PMSM  in  rotor  reference  frame 

Recent  developments  in  microprocessors,  magnetic  materi¬ 
als  and  semiconductors  technologies  have  offered  an  excellent 
opportunity  to  use  ac  motors  in  high  performance  drive  systems. 
PMSM  became  at  the  top  of  ac  motors  in  the  medium  range  of 
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H  =  (Lq  +  /q)/q.  (4) 

Ld  and  Lq  are  the  forward  and  quadratic  axes  cyclic  inductances 
and  l'd  and  l'q  the  forward  and  quadratic  axes  leakage  inductances. 

Expression  of  the  motor  electromagnetic  torque  (Tm)  can  thus 
be  obtained: 


Tm  —  Plq 


^((Ld  +  l'd)  -  (lq  +  l'q))id  - 


(5) 


With  p  the  number  of  poles  and  Vt  the  magnetisation  flux  per 
pole. 

Then  the  variables  such  as  t/t,  e  (electro-motive  force),  Ld,  Tq, 
m  (mutual  between  two  phases  of  machine)  and  the  torque  can 
be  expressed  with  motor  structural  parameters  and  a  completely 
adjustable  model  of  the  motor  is  obtained. 

In  our  case,  1 1  parameters  defining  the  PMSM  structure  are 
considered.  These  parameters  are: 

"  ReCs  (stator  yoke  external  ray)  =  61.6  mm 
hcs (stator  yoke  height)  =  5.9  mm 
h  d(cog  height)  =17.7  mm 
/d(cog  width)  =  5  mm 
ent(air-gap)  =  1  mm 
ha  i  (magnet  height)  =  3  mm 
<  hcr (rotor  yoke  height)  =  6.1  mm 
Lfer(machine  length)  =  30  mm 
p( pole  pair  number)  =  2 
dfii(wire  diameter)  =  1.94  mm 
ns  (windings  number  of  turns)  =  40 
Rs (stator  coil  resistance)  =  0.17  Q 
Ls (stator  coil  inductance)  =  6.83.10-4H 
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Numerical  values  are  the  ones  corresponding  to  the  starting 
values  that  have  been  used  next  in  the  design  process. 

2.2.  Analytical  model  of  compressor 

Rotary  vane  compressor  model  is  constituted  of  four  inputs 
(rotational  speed,  inlet  air  temperature,  and  inlet  and  outlet  air 
pressures)  and  four  outputs  (air  mass  flow,  air  outlet  temperature, 
compressor  torque  and  mechanical  power).  Parameters  used  for 
modelling  are  the  compressed  volume  per  revolution  ( Vcomp/rev) 
obtained  thanks  to  the  compressor  sizes,  air  heat  capacity  (cp) 
and  air  mass  constant  (ra ir). 

The  expression  of  air  mass  flow  ( qm )  is  depending  on  rota¬ 
tional  speed  ( co )  of  compressor  according  to 

Qm  =  t  QV  Vcomp/rev ^Pair  (6) 

Z7T 

with  co:  rated  rotational  speed  (rpm),  pair:  air  inlet  density 
(kg  m-3),  qm :  mass  flow  (kg  s-1),  VComp/rev-  volume  compressed 
by  revolution  (m3rev_1),  qv:  volumetric  efficiency  of  rotary 
vane  compressor. 

Other  important  relation  for  compressor  modelling  is  the 
expression  of  requested  power  linked  to  air  flow  and  pressure 
ratio.  This  relation  is  obtained  from  the  temporal  derivative  of 
the  work  capacity  needed  to  compress  the  air  value  from  the 
inlet  to  the  outlet  pressure. 

1 

P c  —  ^m^pTincomp 

ad  77m 

with  pin  and  pout:  inlet  air  pressure  (Pa)  and  outlet  air  pressure 
(Pa),  7incomp:  inlet  air  temperature  (K),  cy  air  heat  capacity,  r)^: 
adiabatic  efficiency  of  rotary  vane  compressor,  r]m:  mechanic 
efficiency  of  rotary  vane  compressor,  y:  poly  trophic  coefficient 
which  depends  on  compressor  load  conditions  (friction  and  heat 
transfers). 

Finally,  compressor  outlet  pressure  is  defined  by  pressure 
drops  inside  FC  stack  itself.  Experimentally,  it  can  be  verified 
that  the  pressure  drops  are  directly  linked  to  gases  circula¬ 
tion  speed  by  a  coefficient  value  69.3  [2].  Thus  the  following 
empirical  evolution  is  considered  here  (for  the  chosen  FC 
stack): 

Pout  =  Pin  +  69.3  co  (8) 

2.3.  Design  a  fuzzy  logic  controller 

Fuzzy  systems  are  being  used  successfully  in  an  increasing 
number  of  application  areas,  they  use  linguistic  rules  to  describe 
systems.  Fuzzy  logic  provides  a  general  concept  for  description 
and  measurement.  Most  fuzzy  logic  systems  encode  human 
reasoning  into  a  program  to  make  decisions  or  control  a  system. 
Fuzzy  logic  comprises  fuzzy  sets,  which  are  a  way  of  repre¬ 
senting  non- statistical  uncertainty  and  approximate  reasoning, 
which  includes  the  operations,  used  to  make  inferences  in  fuzzy 
logic. 


p  (error  or  error  variation) 


Fig.  2.  Input  and  output  membership  function. 


2.3.1.  The  fuzzy  logic  controller 

The  fuzzy  control  comprises  three  steps:  fuzzification,  infer¬ 
ence,  defuzzification  [8].  In  the  first  step,  the  FLC  receives 
normalized  physical  size  as  input.  The  input  is  then  developed 
in  terms  of  a  set  of  N  membership  functions  (MF)  which  cover 
the  range  of  interest  of  the  corresponding  indicator.  Each  MF  is 
weighted  according  to  its  degree  of  membership  m/  in  the  devel¬ 
opment.  Typically,  for  real-time  control  triangular  MFs  are  used 
for  their  simplicity  (Fig.  2). 

In  the  second  step,  the  inference  operation  represents  the  if- 
then  implication.  The  goal  is  to  interconnect  input  membership 
and  output  membership  of  the  controller.  In  the  last  step,  defuzzi¬ 
fication,  the  results  of  the  inference  step  are  combined  into  a 
single  value  for  each  controller  output. 

In  the  present  case,  for  implementation  of  internal  air  flow 
regulator,  Takagi-Sugeno  fuzzy  logic  controllers  have  been  used 
because  they  allow  minimizing  sampling  time.  In  order  to  reg¬ 
ulate  and  supply  the  required  air  quantity  with  an  acceptable 
response  time  to  the  FC  stacks,  a  non-linear  fuzzy  PD  + 1  reg¬ 
ulator  (Fig.  3)  has  been  chosen  because  of  its  high  dynamic, 
precision  and  robustness  [15]. 

2.3.2.  FLC  considered  optimization  parameters 

In  our  case,  the  following  tuning  parameters  of  the  controller 
are  considered  like  most  significant  (Fig.  3): 

•  Three  parameters  corresponding  to  the  scaling  factors, 

•  Six  parameters  corresponding  to  the  starting  points  and  apex 
points  of  the  membership  functions, 

•  One  parameter  corresponding  to  the  integral  corrective  coef¬ 
ficient. 

Those  are: 

'  PSe  (positive  small  membership  value  on  input  entry  of  controller)  =  0.254 
PVSe (positive  very  small  membership  value  on  error  input)  =  0.033 
PSde (positive  small  membership  value  on  the  error  variation  input)  =  0.7 
PVSde (positive  very  small  membership  value  on  the  error  variation)  =  0.21 
NSs (negative  small  membership  value  on  the  output)  =  —  0.8 
J  NVSs (negative  very  small  membership  value  on  the  output)  =  —0.615 
em(scaling  factor  on  the  error  input)  =  10  g  s_1 

1 

dem(scaling  factor  on  the  error  variation  input)  =  — (r  +  0AT)Teem 

x.T 

2.07 

gm(scaling  factor  on  the  output)  =  - (r  +  0.4T)em 

KT 

1.6 

^•(coefficient  of  the  integral  correction)  =  - 


Poxxi  A 
Pm  ) 


(y-i/y) 


-  1 
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Fig.  3.  Scheme  of  fuzzy  PD  + 1  controller  principle. 

with  r:  time-constant,  T:  delay,  K:  static  gain,  TQ:  sampling 
period. 


3.  Optimization 

3.1.  Considered  operating  cycle  and  objective  function 


The  complete  air  supply  circuit  system  is  tested  on  a  step 
set  point  and  on  a  vehicle  driving  cycle.  The  first  one  presents 
an  amplitude  of  10  gs-1  during  a  138-ms  cycle  and  the  sec¬ 
ond  one  is  a  vehicle  normalized  mission  profile  (Fig.  4)  that  is 
deduced  from  the  vehicle  simulator  program,  Advisor®  (profile 
“CYC_UDDS:  Urban  Dynamometer  Driving  Schedule”),  used 
in  embedded  FC  systems  experimentation  [16]. 

In  order  to  optimize  this  system  for  these  set  points  the  objec¬ 
tive  function  has  to  be  defined.  For  example,  objective  function 
(/)  is  here,  for  step  set  point  the  following: 


/  = 


rt  ec 

rrt  •  2  rec  •  2 


(9) 


with  rt  the  rising  time,  rrt  the  reference  system  rising  time,  ec 
the  energy  consumption  and  rec  the  reference  system  energy 
consumption. 

That  allows  obtaining  a  50%  influence  of  the  energy  crite¬ 
rion  and  a  50%  influence  of  the  rising  time  on  the  optimization, 
leading  to  both  an  economic  and  dynamic  system. 


3.2.  The  reference  system 


Classically,  optimization  methodologies  applied  to  this  kind 
of  system  are  sequential,  firstly  optimizing  the  actuator  alone, 


x  icr3 


Fig.  4.  Fuel  cell  air  mass  flow  (kgs  ])  consumption  profile  deduced  from 
Advisor®  (CYC_UDDS)  on  60s. 


then,  defining  the  control  for  the  actuator.  Nevertheless,  a  lim¬ 
itation  of  this  approach  is  that  the  capacities  of  the  control  are 
then  limited  by  the  actuator  structure  itself.  The  proposed  new 
approach  carries  out  system  global  optimization  in  calculating, 
jointly,  both  the  optimal  actuator  structural  parameters  and  the 
optimal  control  law  parameters. 

Thanks  to  the  values  of  the  reference  system  (given  in  pre¬ 
vious  section),  the  general  shape  of  the  motor  (Fig.  5)  and  the 
control  surface  (Fig.  6),  that  reflects  the  control  strategy  imposed 
on  the  controller,  can  be  obtained.  This  control  surface  gives  a  3D 
representation  of  the  control  law  realized  by  the  FLC.  The  two 
first  dimensions  are  the  two  normalized  inputs  of  the  FLC,  the 
last  dimension  gives  the  normalized  control  output  of  the  FLC. 
The  reachable  power  greatly  depends  on  the  air  flow  within  the 
stack.  If  the  air  supply  circuit  cannot  respond  quickly  enough 
to  the  requested  air  flow,  the  FC  life  time  is  restricted.  Thus,  a 
chosen  air  flow  dynamic  (rise  time:  138  ms)  is  here  supposed  to 
be  requested  by  the  FC.  This  rise  time  value  has  been  obtained 
considering  the  oxygen  volume  (0.69  g)  that  is  present  inside  the 
cathode  compartment  and  the  time  during  which  this  oxygen  can 
supply  the  fuel  cell  in  case  of  a  nominal  current  set  point  on  our 
real  experimental  5  kW  FC  stack  [17].  With  this  system,  the  air 
flow  responses,  obtained  with  the  initial  parameters  settings,  are 
plotted  in  Fig.  7  for  step  set  point  and  Fig.  8  for  a  mission  profile. 

According  to  Fig.  7,  the  rising  time  is  125  ms  and  no  steady 
static  error  is  observed  (for  the  reference  system).  The  require¬ 
ment  of  a  rising  time  lower  than  138  ms  is  respected.  The  system 
consumes  764  J  for  this  cycle. 


0.5 
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-0.5 
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Fig.  6.  The  reference  system  control  law  strategy. 
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Fig.  7.  The  reference  system  response  of  air  mass  flow  to  step  set  point. 
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Fig.  9.  Principle  of  the  displacement  of  a  particle  in  a  two-dimensional  space. 


Exploration  domain  for  each  variable  with  Recs,  hcs ,  hd.  Id,  ent, 
ha i,  hcv,  Lfer,  dm  in  millimeter. 


Fig.  8  gives  the  response  of  the  system  for  dynamical  solici¬ 
tation.  We  can  see  that  smearing  errors  is  almost  null  (maximum 
transient  error  is  about  8.17  x  10-4gs-1).  The  system  con¬ 
sumes  45,129  J  for  this  cycle  of  60  s. 

Now  let  us  define  the  exploration  domain  for  each  optimiza¬ 
tion  parameter. 


RecsC  [0.01,  0.1] 
hcse  [0.003,  0.01] 
hde  [0.005,  0.03] 
lde  [0.003,  0.01] 
ent  g  [0.0005,  0.005] 
/lai  £  [0.003,  0.01] 
hCY  €  [0.003,  0.015] 
Lfere  [0.01,  0.1] 

pel  1,6] 

dm  g  [0.0002,  0.005] 
I  ns  g  [0,  100] 


^  PSe  G  [0,  1] 
PVSe  G  [0,  1] 
PSde  €  [0,  1] 
PVSde  €  [0,  1] 

<  NSs  G  [—1,  0] 
NVS,g[-1,0] 
de  g  [0,  10] 
^€[0,  10] 

I  k  G  [0,  10] 


(10) 


3.3.  Particle  swarm  optimization 
3.3.1.  Laws  of  the  PSO 

PSO  belongs  to  the  broad  class  of  stochastic  optimization 
algorithms.  The  ideas  that  underlie  PSO  are  not  inspired  by  the 
evolutionary  mechanisms  encountered  in  natural  selection,  but 
rather  by  the  social  behaviour  of  flocking  organisms,  such  as 
swarms  of  honeybees  and  fish  shoals.  It  has  been  observed  that 
the  behaviour  of  the  individuals  that  belong  to  a  flock  adheres 
to  fundamental  rules  like  nearest-neighbour  velocity  matching 
and  acceleration  by  distance  [18,19].  PSO  is  a  population-based 
algorithm  that  uses  a  population  of  individuals  to  probe  promis¬ 
ing  regions  of  the  search  space.  In  this  context,  the  population  is 
called  a  swarm  and  the  individuals  are  called  particles.  Each  par¬ 
ticle  moves  with  an  adaptable  velocity  within  the  search  space 
(Fig.  9),  and  keep  in  its  memory  the  best  position  it  ever  encoun¬ 
tered.  In  the  global  variant  of  PSO,  the  best  position  ever  reached 
by  all  individuals  of  the  swarm  is  communicated  to  the  whole 
particles.  In  the  local  variant,  each  particle  is  assigned  to  a  neigh¬ 
bourhood  consisting  of  a  specified  number  of  particles.  In  this 


Fig.  8.  The  reference  system  response  of  air  mass  flow  to  step  dynamic  solicitation. 
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case,  the  best  position  ever  reached  by  the  particles  that  com¬ 
prise  the  neighbourhood  is  communicated  among  them  [20]. 
This  paper  considers  the  global  variant  of  PSO  only. 

Let  us  assume  a  D-dimensional  search  space  S,  Sc  IRD, 
and  a  swarm  consisting  of  N  particles.  The  ith  particle  is  a 
D-dimensional  vector  has  the  following  coordinates: 

=  (Xi\  i  A’ilt  •  •  •  ?  Xijf)  £  S 

The  velocity  of  this  particle  is  also  a  D-dimensional  vector  such 
as 

Vi  =  ( vn ,  va, .  ■  ■  viD)T  e  S 

The  best  previous  position  encountered  by  the  ith  particle  is  a 
point  in  S ,  denoted  as 

Pi  =  (pn,  pa,  ■  ■  ■  PiD)T 

If  g  is  the  index  of  the  particle  that  reached  the  best  previous 
position  among  all  the  individuals  of  the  swarm,  and  t  is  the  iter¬ 
ation  counter.  Then,  according  to  the  version  of  PSO,  a  parameter 
called  the  constriction  factor  is  defined  by  Clerc  [21] .  The  swarm 
is  thus  manipulated  according  to  the  following  equations: 

Vi(t  +  1)  =  -  Xi(t)) 

+  c2r2(Pg(t)-Xi(t ))]  (11) 


Xi(t  +  1)  =  Xi(t)  +  Vi(t  +  1)  (12) 


employed  a  Pareto  ranking  scheme  to  handle  constraints,  which 
is  a  concept  of  multi-objective  optimization. 

Here,  a  simple  but  efficient  method  is  introduced  to  solve 
constrained  optimization  problems.  The  preserving  feasibility 
strategy  is  implemented  to  deal  with  constraints.  Two  modifica¬ 
tions  of  the  PSO  algorithms  have  been  proposed: 

(1)  when  updating  the  memories  pBest  (coordinates  of  best  fit¬ 
ness  value  of  each  particles)  and  gBest  (coordinates  of  best 
fitness  value  of  all  particles),  each  particle  only  keep  feasi¬ 
ble  solutions  in  its  memory  and  leave  aside  other  solutions 
that  do  not  respect  imposed  constraints, 

(2)  during  the  initialization  process,  one  particle,  at  least,  is 
initialized  started  with  feasible  solution. 

Compared  to  other  constraint  handling  techniques,  this 
approach  has  the  following  advantages: 

(1)  It  is  quite  simple.  There  is  no  pre-processing  to  take  into 
account  the  constraints  and  there  is  neither  complicated 
manipulation.  Fitness  function  and  constraints  are  handled 
separately,  thus  there  are  no  limitations  concerning  the  con¬ 
straints. 

(2)  It  is  faster.  The  only  part  of  the  algorithm  dealing  with  con¬ 
straints  is  to  check  if  a  solution  satisfies  all  the  constraints. 
This  will  reduce  the  computation  time  when  handling  mul¬ 
tiple  or  complicated  constraints. 

Examples  of  physical  constraints  are  given  as  following: 


•  External  radius  of  the  rotor  yoke  <  external  radius  of  the  magnets, 

•  Diameter  of  winding  wire  <  spaces  passage  between  two  isthmuses  of  tooth, 

•  PSe  >  PVSe,  1 

•  PSde  >  PVSde,  J  Fuzzy  loSlc  parameters 

•  35%  >rate  of  filling  of  the  notches  ... 


where  i=  1,  2,  . . .,  N,  x  is  the  constriction  factor,  c\  and  c2 
(with  experimental  value  of  2.5)  denote  the  cognitive  and  social 
parameters,  respectively,  and  r\,  r2  are  random  numbers  uni¬ 
formly  distributed  in  the  interval  [0,1]. 

The  value  of  the  constriction  factor  is  typically  obtained 
through  the  formula  [21]: 

2k 

X  =  -  ,  for  0  >  4 

0  —  2  +  sj 02  —  40  with  k  =  1 

X  =  yfic,  for  0  <  4 


3.4.1.  Implementation  of  PSO 

The  implementation  of  PSO  program  is  easy  and  only  takes  a 
few  lines.  For  example,  in  our  case,  the  algorithm  programming 
on  Matlab®  uses  25  lines.  More  precisely  the  different  steps  of 
the  PSO  program  are  described  in  Fig.  10. 


and  0  =  ciri  +  c2r2 


>(13) 


3.4.  Constraint  handling 

In  the  literature,  several  studies,  which  proposed  to  extend 
PSO  to  constrained  optimization  problems,  are  reported,  and 
different  constraint  handling  techniques  were  used. 

Parsopoulus  [22]  converted  the  constrained  optimization 
problem  into  a  non-constrained  optimization  problem  by  adopt¬ 
ing  a  non-stationary  multi-stage  assignment  penalty  function 
and  then  applying  PSO  to  the  converted  problems.  Ray  [23] 


3.5.  Optimization  process 

The  control  concept,  using  field  oriented  control  and  imple¬ 
mentation  of  fuzzy  logic  controller  on  the  motor-compressor 
group  is  shown  in  Fig.  11.  Firstly,  for  this  model,  the  PSO 
algorithm  generates  motor  construction  parameters  and  con¬ 
trol  law  parameters.  In  the  second  step,  the  model  is  simulated. 
The  third  step  consists  in  simulation  result  evaluation  in  order 
to  calculate  the  value  of  optimization  criterion  /.  Finally, 


A.  Sari  et  al.  /  Journal  of  Power  Sources  179  (2008)  121-131 


127 


Fig.  10.  Flow  chart  of  optimization  process  by  PSO. 


a  new  set  of  optimization  parameters  is  evaluated  (cf.  step 

1). 

In  our  case,  all  optimizations  are  based  on  401  iterations  and 
20  particles  in  order  to  ensure  equity  in  comparing  the  methods. 
According  to  numerous  papers,  the  choice  of  20  particles  is 
suitable  for  a  wide  number  of  optimization  processes.  Stochastic 
optimization  methods  ensure  finding  the  global  optimum  with  an 
infinite  time  of  optimization  what  is  not  really  conceivable.  The 
choice  of  401  iterations  is  based  on  our  own  simulation  know¬ 
how  and  on  the  constraint  to  have  a  reasonable  computational 
time. 

Let  us  have  a  look  at  the  optimization  progress  (cf. 
Figs.  12  and  13).  On  Fig.  12,  the  progress  (during  iterations) 
of  the  state  of  particles  is  provided.  The  particles  respecting  all 
the  imposed  constraints  are  entered  as  good  particles.  We  can 
observe  the  increase  in  the  number  of  good  particles  with  the 


fitness 


Fig.  12.  Evolution  of  good  particles  percentage  during  optimization. 


increase  in  the  iteration  count.  Reaching  a  certain  iteration  count 
we  see  a  sudden  fall  of  the  number  of  good  particles  which  can 
be  explained  by  the  process  of  dispatching,  i.e.  when  90%  of  the 
particles  are  in  an  hypersphere  dimension  R  (another  variable), 
they  are  dispersed  randomly  in  search  space.  That  is  a  comple¬ 
mentary  function  making  it  possible  to  decrease  the  probabilities 
of  being  blocked  in  a  local  minimum. 

Fig.  13  provides  the  evolution  of  the  objective  function  which 
begins  with  a  value  of  0.807  that  corresponds  to  the  objective 
function  obtained  for  the  reference  system  and  stops  with  a  value 
of  0.31  that  is  the  value  of  the  objective  function  obtained  on 
the  system  defined  by  simultaneous  optimization  on  an  indi- 
cial  solicitation.  The  objective  function  evolution  becomes  very 
slow  after  250  iterations.  Furthermore,  particles  are  three  times 
redistributed  between  iteration  250  and  iteration  401.  It  is  thus 
possible  to  say  it  is  at  least  a  local  optimum. 

4.  Results 

4.1.  Sequential  optimization 

Firstly,  a  sequential  optimization  method  for  the  motor- 
compressor  and  its  control  has  been  carried  out.  The 
optimization  process  starts  with  PMSM  optimization  based  on 
an  efficiency  criterion.  PSO  leads  to  a  PM  motor  with  an  effi¬ 
ciency  of  95%  (Fig.  14).  Then  using  PSO,  optimal  control  law 
parameters  are  computed  (considering  the  two  types  of  set  points 
that  have  been  considered).  For  instance,  considering  a  step  set 


Ambient  air 


Fig.  11.  System  model. 
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best  value  profile 


Fig.  13.  Objective  function  evolution  during  optimization. 


Fig.  14.  PMSM  structure  (for  indicial  set  point)  for  the  system  optimized  by 
sequential  method. 

point,  the  control  law  obtained  by  PSO  enable  us  to  build  the 
control  surface  presented  in  Fig.  15. 

This  motor-compressor  group  system  provides  a  response  in 
mass  flow  with  a  rising  time  of  132  ms,  without  steady-state 
error  and  with  a  consumption  of  459  J  for  this  cycle  (Fig.  16). 

Now,  for  the  dynamic  solicitation  the  control  law  obtained  by 
PSO  enables  to  build  the  control  surface  of  Fig.  17.  These  control 
law  and  motor  design  provide  a  response  in  mass  flow  with  very 
weak  smearing  errors  corresponding  to  a  maximal  interval  of 
0.689  and  a  consumption  of  9718  J  on  this  60  s  cycle  (Fig.  18). 


Fig.  15.  Control  law  surface  (for  indicial  set  point)  for  the  system  optimized  by 
sequential  method. 


Fig.  17.  Control  surface  obtained  by  sequential  optimization. 


4.2.  Coupled  optimization 

Secondly,  simultaneous  optimization  method  for  the  motor- 
compressor  and  its  control  has  been  performed.  Using  this 
method,  PSO  carries  out  an  iterative  research  in  order  to  obtain 
the  21  parameters  of  the  engine  and  the  fuzzy  controller  param¬ 
eters  that  will  minimize  the  objective  function.  The  objective 
function  is  defined  to  achieve  the  goals  of  desired  consumption 
and  performance. 


Fig.  18.  Response  in  air  mass  flow  for  the  dynamic  set  point. 
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Fig.  20.  Control  law  surface  for  the  system  optimized  simultaneously. 


For  the  step  set  point,  the  parameters  obtained  by  PSO  gen¬ 
erate  the  motor  geometry  and  the  control  surface  shown  in 
Figs.  19  and  20. 

In  this  case,  a  rising  time  (for  air  mass  flow)  of  25  s  and  a 
steady-state  error  are  obtained  with  a  light  over- shooting,  and  a 
consumption  of  568  J.  We  can  see  the  response  in  Fig.  21. 

In  short,  we  obtain,  with  a  simultaneous  optimization,  a  sys¬ 
tem  much  more  powerful  than  the  reference  system  or  the  system 
optimized  sequentially  and  a  system  that  has  power  consump¬ 
tion  much  lower  than  the  reference  one  but  a  little  energy-hungry 
than  the  sequentially  optimized  system. 


Fig.  21.  Response  in  air  mass  flow  step  set  point. 


Fig.  22.  PMSM  for  the  system  optimized  by  simultaneous  optimization. 

The  greediness  of  the  system  optimized  simultaneously 
comes  from  the  ratio  between  the  energy  consumption  and  the 
performance  to  form  the  objective  function.  Thus  by  modifying 
this  ratio,  different  systems  could  have  been  obtained  (less  or 
more  dynamic,  less  or  more  powerful). 

For  the  dynamic  set  point,  the  parameters  given  by  PSO  gen¬ 
erate  the  engine  profile  and  the  control  surface  of  Figs.  22  and  23. 

These  control  laws  and  engine  give  us  a  response  in  mass  flow 
with  a  very  small  smearing  errors  (maximal  interval  of  0.4788  g) 
and  a  consumption  of  9133  J  for  this  cycle  of  60  s.  The  response 
is  shown  in  Fig.  24. 

Thus,  with  a  simultaneous  optimization,  we  obtain  a  system 
much  more  economic  than  the  sequentially  optimized  system 
and  on  enormously  more  economic  than  the  reference  system. 
The  counterpart  is  that  the  motor  obtained  by  sequential  opti¬ 
mization  is  double  sized  versus  the  reference  one  and  that  the 
machine  obtained  by  simultaneous  optimization  is  still  a  little 
bigger  than  the  last  one  (Table  1).  Ultimately,  the  previous  results 
prove  the  effectiveness  of  the  simultaneous  optimization  method 
using  the  particular  swarms. 

4.3.  Coupled  optimization  with  mass/energy  criterion 

Previous  optimizations  allow  seeing  PSO  sizing  capacity  and 
improvement  ability  by  system  approach.  But  to  obtain  a  real¬ 
istic  and  useful  structure  of  PMSM,  it  is  necessary  to  adapt  the 
optimization  criterion.  When  we  optimize  only  on  an  energetic 


Fig.  23.  Control  law  surface  for  the  system  optimized  by  simultaneous  opti¬ 
mization. 
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Fig.  24.  Response  in  air  mass  flow  at  the  dynamic  set  point. 


Table  1 


Summary  table  of  the  optimized  system  characteristics 


Characteristic  size 

Initial  system 

System  optimized  sequentially 
for  a  dynamic  solicitation  and  an 
energy  criterion 

System  optimized  simultaneously 
for  a  dynamic  solicitation  and  an 
energy  criterion 

System  optimized  simultaneously 
for  a  dynamic  solicitation  and  a 
mass/energy  criterion 

External  radius 

61.6  mm 

79.1  mm 

70.4  mm 

75.7  mm 

Motor  length 

30  mm 

49.2  mm 

73.9  mm 

16.05  mm 

Pole  pair 

2 

2 

2 

6 

Energy  consumption 

45,129  J 

9,718  J 

9,133  J 

10,886  J 

Motor  mass 

3.39  kg 

6.94  kg 

8.72  kg 

1.07  kg 

criterion,  the  algorithm  tends  to  increase  motor  length  to  mini¬ 
mize  energy  consumption.  But  it  is  not  always  acceptable.  It  is, 
thus,  more  interesting  to  take  an  optimization  criterion  which 
is  a  compromise  between  energetic  consumption  and  mass  of 
PMSM.  Optimization  on  401  iterations  with  20  particles  in  a 
unique  swarm  with  a  mass/energy  criterion  allows  obtaining  one 
new  system.  Fig.  25  shows  PMSM  structure  and  Fig.  26  shows 
fuzzy  logic  control  law  surface.  This  new  system  has  an  answer 
(Fig.  27)  for  the  dynamic  solicitation  with  a  maximal  smearing 
error  of  4.34e-5  g  and  a  consumption  of  10,886  J  on  the  cycle 
of  60s.  This  system  has  more  important  energy  consumption 
that  both  systems  previously  optimized  but  the  energy  motor  is 


Fig.  25.  PMSM  for  the  system  optimized  by  simultaneous  optimization  with  a 
mass/energy  criterion. 


Fig.  26.  Control  law  surface  for  the  system  optimized  by  simultaneous  opti¬ 
mization  with  a  mass/energy  criterion. 

x  10'3 


Fig.  27.  Response  in  air  mass  flow  for  the  dynamic  set  point  for  the  system 
optimized  by  simultaneous  optimization  with  a  mass/energy  criterion. 
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smaller.  Table  1  compares  obtained  system  characteristics.  It  is 
sharply  visible  that  for  about  20%  more  important  energy  con¬ 
sumption,  motor  mass  is  divided  by  8.  For  an  embedded  system, 
this  mass  gain  is  very  consequent. 

5.  Conclusion 

In  this  paper,  an  innovating  optimization  approach  for  the 
design  of  a  FC  motor-compressor  group  shows  its  effective¬ 
ness  based  on  the  results  obtained  for  various  design  tests.  This 
method,  known  as  simultaneous  design,  uses  a  system  approach 
in  order  to  design  the  whole  system.  This  method  combined 
with  a  heuristics  such  as  PSO,  shows  its  effectiveness  in  terms 
of  implementation  simplicity  and  required  computing  time.  The 
different  results  proved  that  simultaneous  optimization  gives 
better  results  on  this  system  than  the  sequential  approach.  How¬ 
ever  this  approach  must  be  validated  on  experimentation  in  real 
time.  Moreover,  it  is  presumable  to  still  improve  our  results  by 
including  the  parameters  of  the  compression  head  model  in  the 
simultaneous  optimization.  Indeed,  it  will  give  a  greater  degree 
of  freedom  to  the  optimization  program  in  order  to  achieve 
our  goal  to  minimize  the  power  consumption  of  the  motor- 
compressor  on  a  given  mission  cycle.  Improvement  could  also  be 
made  into  the  analytical  model  of  PMSM  by  taking  into  account 
excess  losses  such  as  iron  or  mechanical  losses.  This  work  is  on 
progress  in  our  laboratory. 
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